Thermodynamics of the quantum easy-plane antiferromagnet on the triangular lattice 
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The classical XXZ triangular-lattice antiferromagnet (TAF) shows both an Ising and a BKT tran- 
sition, related to the chirality and the in-plane spin components, respectively. In this paper the 
quantum effects on the thermodynamic quantities are evaluated by means of the pure-quantum 
self-consistent harmonic approximation (PQSCHA), that allows one to deal with any spin value 
through classical MC simulations. We report the internal energy, the specific heat, and the in-plane 
correlation length of the quantum XXO TAF, for S = 1/2, 1, 5/2. The quantum transition temper- 
atures turn out to be smaller the smaller the spin, and agree with the few available theoretical and 
numerical estimates. 



A renewed interest has recently focused on triangu- 
lar antiferromagnets (TAF) ^. Indeed, they turned 
out to describe the magnetic behavior of several real 
compounds as, for example, the stacked antiferromag- 
net NaTi02 [||, the organic superconductors of the fam- 
ily K-(BEDT-TTF)2X || and the K/Si(lll):B inter- 
face I]. 

In this paper we investigate the thermodynamic prop- 
erties of the quantum XXZ Heisenberg antiferromagnet 
on the triangular lattice, defined by the following Hamil- 
tonian 
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where J is the positive (antiferromagnetic) exchange con- 
stant, and (5?, Sy , S?) are the spin operators sitting on 
the sites i of a triangular lattice. They satisfy SU{2) 
commutation relations and belong to the spin-S" repre- 
sentation |S?| = S{S -\- 1). The interaction is restricted 
to nearest-neighbors and d runs over their relative dis- 
placements. The planar character of the system is due to 
the presence of the anisotropy A € [0, 1), energetically fa- 
voring configurations with the spins lying in the xy plane 
(easy-plane). For A = the spin components on the z 
axis do not appear in the Hamiltonian and the model is 
known as XXO or quantum XY. 

The minimum energy configuration of the classical 
counterpart of the Hamiltonian (nl) for every value of 
A £ [0, 1] consists of coplanar spins forming ±27r/3 angles 
between nearest-neighbors and this leads to a a/3x \/3 pe- 
riodic Neel state. In contrast to the isotropic case, where 
the plane in which the 27r/3 structure lies can take any 
direction in the spin space, in the XXZ model such struc- 
ture must take place in the easy-plane. As a result, in the 



planar TAF the frustration causes an additional discrete 
two-fold degeneracy of the classical ground state, which 
is due to the chirality (or helicity), defined as the sign of 
rotation of the spins along the sides of each elementary 
triangle. The resulting degeneracy corresponds to the 
group S0{2) X Z2. As the Mermin and Wagner theorem 
only states that the sublattice magnetization must van- 
ish at any non zero temperature, long-range order can 
occur as far as the chirality is concerned, and an Ising- 
like phase transition is indeed observed g , in addition to 
the usual Bcrezinskii-Kosterlitz-Thouless (BKT) critical 
behavior associated to the rotation symmetry in the xy 
plane. 

In the quantum case the situation is far less clear. 
In fact, unlike the antiferromagnet on the square lat- 
tice where there is a general consensus about the or- 
dered nature of the ground state even for S = 1/2, in 
the frustrated cases the lack of exact analytical results 
is accompanied by difficulties in applying stochastic nu- 
merical methods, as their reliability is strongly limited by 
the well-known sign problem. Indeed only very recently 
a systematic size-scaling of the order parameter and of 
the spin gap has been performed using a new Quantum 
Monte Carlo technique q], confirming the existence of 
Neel long-range order in the ground state as also sug- 
gested by the symmetry properties of the first excited 
states, evidenced by Bernu et al. Q]. 

An even less clear situation is that concerning the fi- 
nite temperature behavior. In fact an early numerical 
work Ig], limited to lattice sizes up to 27 sites, indicated 
for the S = 1/2 XXO model a phase diagram similar 
to the classical one, in contrast with the high tempera- 
ture expansion produced by Fujiki and Belts [^ where 
no evidence for a phase transition was found. For the 



XXZ Haniiltonian, Momoi and Suzuki |TQ], applying an 
effective field theory, conjectured that the chiral phase 
transition should persist for every value of A € [0, 1), as 
in the classical case, and were able to estimate the tran- 
sition temperature for A = 0, obtaining a value very close 
to that found in Ref. M. Recently Suzuki and Matsub- 
ara [0| using a quantum transfer Monte Carlo method 
to study clusters up to 24 sites, have claimed instead the 
absence of the chiral order at any finite temperature for 
A>0.6. 

In this context where, at least up to now, quan- 
tum numerical methods cannot give satisfactory answers, 
the pure-quantum self- consistent harmonic approxima- 
tion (PQSCHA) ||lj] can provide an effective instrument 
to investigate the thermodynamics of quantum spin sys- 
tems, as far as their ground state is ordered. The method 
is based on the path-integral formulation of quantum sta- 
tistical mechanics, and has been successfully applied re- 
cently to a variety of unfru strate d spin models, both one- 
||13| and two-dimensional 
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By the PQSCHA the evaluation of thermal averages 
in the quantum model can be reduced to the calcula- 
tion of classical-like averages over a Boltzmann distri- 
bution defined by an effective Haniiltonian, which con- 
tains the contribution of the pure-quantum part of the 
fluctuations (approximated within a self consistent har- 
monic scheme) in its renormalized interaction parame- 
ters, which are temperature and spin dependent. As a 
result one can get accurate results on the quantum spin 
system using classical computational methods, like the 
transfer-matrix in the one-dimensional case and classical 
Monte Carlo simulations in the two-dimensional one. 

The first step of the derivation of the effective Hamil- 
tonian for the easy-plane TAF, is to apply the unitary 
transformation which defines a spatially varying coor- 
dinate system pointing along the local Neel direction, 
namely 



U = exp 
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where B and C labels two of the three sublattices. Unlike 
in the bipartite lattices where the corresponding trans- 
formation maps the antiferromagnet into a model with 
an in-plane ferromagnetic exchange and an antiferromag- 
netic coupling along the z axis, thus allowing the demon- 
stration of the Lieb and Mattis theorem ||l6[ and com- 
putability with standard quantum Monte Carlo meth- 
ods, in the triangular case the transformed Hamiltonian 
shows an extra current-like term [|^ which contains the 
effects of the frustration, whose form is quite similar to 
the chiral order parameter Q, i.e., the physical quantity 
undergoing the order-disorder phase transition present in 
the classical model for every value of A < 1. 

From now on the derivation follows the same lines al- 
ready described in Refs. 0,15[. A point worth being 



recalled is the use of the Villain transformation ||I^ in 
order to represent the spin operators in terms of canon- 
ically conjugated variables, which is a necessary step 
in the derivation of the effective Hamiltonian. As it 
is well known, this spin-boson transformation preserves 
the commutation rules but neglects the so called kine- 
matic interaction due to the limited spectrum of S^, 
thus giving a better description when the system has 
a good easy-plane character and the spin states with 
large fluctuations of S^ are less relevant to the thermo- 
dynamics. In the square lattice case Pq | such approxi- 
mation scheme turns out to be reliable up to some value 
of Am < 1 (Am = 0.58 in the extreme quantum case 
S = 1/2), when the mapping with the Villain transfor- 
mation breaks down and a different spin-boson transfor- 
mation is needed. However, it provides accurate results 
for the critical temperatures even for A = 0.5; a similar 
behavior is also found in the case of the quantum TAF. 
Finally we remind that Weyl ordering, which is inherent 
to the PQSCHA, naturally leads to define an effective 
classical spin length a.s S — S -\- 1/2 and thus to set the 
natural energy scale e = JS^. Therefore in the following 
we use the reduced temperature, t = ksT j J S^ . 

In the case of the easy-plane TAF the effective Hamil- 
tonian has the form TicS = 'H + G{t), where G{t) is an ad- 
ditive uniform term, formally identical to that obtained 
for the square lattice, which is unessential in the calcula- 
tion of the thermal averages, while: 
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where (sfjsfjSf) are unit vectors, i.e., classical spins. 
Within the PQSCHA quantum effects are embodied in 
the temperature and spin dependence of the renormalized 
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with 
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where 



a^ =-(l-ii?^)e-2^ii (l + 2Aeff7k) , (8) 
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/k = ak&k/(2S'i) , C{x) = cothx — x~^ is the Langevin 
function, 7k = •^"^X^d cos(k-d) , and k is a wavevec- 
tor varying in the first Brillouin zone. D±{S, X,t) and 
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FIG. 1. Temperature and spin dependence of the inter- 
nal energy per spin for A = and (from the top curve) 
S = 1/2,1,3/2,5/2,5 and oo. SoUd circles are classical MC 
data. In the inset the derivative of the effective exchange 
constant is plotted vs. t with the same convention on the 
lines. 

I?ll (5, A, t) represent the pure-quantum square fluctua- 
tions of the out-of-plane and in-plane components of the 
spins respectively. They are decreasing functions of tem- 
perature and spin, vanishing for i — > oo and S —^ oo, i.e., 
when the quantum part of the fluctuations is negligible 
with respect to the classical one. 

From the above equations, we can infer that the 
PQSCHA approach is valid under the condition that sec- 
ond order terms in {D±)'^ can be neglected. One can take 
the criterion that the renormalization effects of quantum 
fluctuations must not reduce much more than, say, 50% 
the effective exchange integral. Such a strong renormal- 
ization only occurs for S = 1/2 and t < 0.2, while for 
higher spin values the PQSCHA is reliable at any tem- 
perature. 

In the XXO model, Aofr = A = and all the informa- 
tion about the quantum system is hence contained in the 
renormalization of the energy scale. In this case the crit- 
ical properties of the quantum system at a temperature 
t are essentially those of its classical counterpart at the 
effective temperature icff = t/jos(t, S), and we have used 
the results of classical Monte Carlo (MC) simulations re- 
cently obtained for lattice sizes up to iV = 120 x 120 g] 
to calculate the corresponding quantum observables. 

Indicating the classical averages with the effective 
Hamiltonian as (. . ■)^ff, the internal energy per spin can 
be calculated as 



e(t,5,A) = M^(H)^^ 
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where 'D± can be expressed as D± , Eq. (o) , with an extra 
factor 7k in the summand. For A = the above equation 
reduces to 




FIG. 2. Temperature and spin dependence of the specific 
heat in the XXO model. Circles are the classical MC data 
obtained from the mean squared fluctuations of the energy 
while the solid line is the numerical derivative of the energy 
curve. 



6(^,5") = Jeff (i) eel (ioff), 
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being eei(t) the internal energy per spin of the corre- 
sponding classical system. In Fig. |l| e{t, S) is plotted for 
various values of the spin in the range of temperatures 
where the PQSCHA is expected to give reliable results. 

The energy curves flatten and increase with decreasing 
S due to the increased quantum fluctuations. As said be- 
fore the 5 = 1/2 curve is reported only in the valid tem- 
perature range. As a matter of fact the extrapolation to 
lowest temperatures gives the self-consistent spin-wave 
ground-state energy. The difference from the most re- 
fined estimates |1^ can be mainly attributed to 1/5'^ 
constant contributions coming from the Villain transfor- 
mation Ml and to the use of the low coupling approxima- 
tion (LCA) ||T^. This term is not significant for S* > 1. 

Consistently with this picture the finite size {N = 
120 X 120) peak of the specific heat (Fig. ||), obtained 
by numerical derivation of the internal energy, moves 
towards lower temperatures and decreases in height as 
S decreases. However, since the quantum renormaliza- 
tions are essentially size- independent, classical scaling 
with size pj is conserved and a logarithmic divergence 
of the specific heat, connected with the Ising-like chiral- 
ity phase-transition, is therefore expected in the ther- 
modynamic limit. By direct derivation of Eq. (|ll| ) it is 
easily seen that, in order for the quantum specific heat to 
vanish in the zero temperature limit, within our approx- 
imation we must have djcs/dt —^ |eci(0)|~^ as t — > a 
condition which is fulfilled for every value of S, as can be 
verified analytically from the explicit expressions of the 
renormalization parameters; this is shown in the inset of 
Fig.0. 

Most papers in literature do mainly concern with the 
chiral order-disorder transition, while the XXZ TAF also 
supports another kind of phase transition. In fact, the 




FIG. 3. Temperature and spin dependence of the magnetic 
correlation length in the XXO model. Circles and solid line 
are classical MC data. 



classical system g displays as well a BKT critical be- 
havior. For this reason we have calculated the mag- 
netic correlation length which governs the decay of the 
in-plane correlation functions in the high-temperature 
phase, whose expression within the PQSCHA reads 
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where i and i+r belong to the same sublattice and G{r) 
is bounded and essentially constant for large r. As a 
consequence the asymptotic behavior of the correlation 
functions in the critical region is the same of the effective 
classical spin system. In particular, for A = 0, the corre- 
lation length can be simply found as £,{t) = Cci(^cff)- The 
result is reported in Fig. y: as expected in a BKT tran- 
sition, it displays a divergence at a temperature isKT 
which decreases with decreasing 5", as a result of en- 
hanced quantum fluctuations. 

Within the PQSCHA, the quantum renormalizations 
cannot modify the critical behavior of the effective clas- 
sical system g^ , so that both the chirality and the BKT 
critical temperatures can be connected to their classical 
counterparts by the self-consistent relation 
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which can be solved numerically. The obtained critical 
temperatures for A = and A = 0.5 arc reported for 
various values of the spin in Table |. In the S = 1/2 case, 
although our theory begins to become unreliable when 
t < 0.2, we notice that for A = the extrapolated value 
for the chiral critical temperature, tc = 0.193, agrees 
remarkably well with those obtained by the size scaling 
on the QMC § data (tc = 0.195(1)) and the effective 
field theory of Ref. @ (ic ^ 0.20). 

L. C. acknowledges S. Sorella for continuous and fruit- 
ful discussions. 



TABLE I. Chiral and BKT critical temperatures for A = 
and A — 0.5 and for some values of the spin length S. The 
classical values are taken from Ref. M. The reported errors 
only represent the statistical uncertainty of the MC data. 



1/2 



3/2 



5/2 



tciS,0) 0.193(2) 0.273(3) 0.319(3) 0.364(4) 0.396(5) 0.412(5) 

tBKT(S, 0) 0.1875(5) 0.265(1) 0.310(1) 0.352(1) 0.386(1) 0.402(2) 

ic(S, 0.5) 0.185(4) 0.267(4) 0.312(5) 0.355(5) 0.385(5) 0.400(5) 

tBKT(S, 0.5) 0.180(1) 0.260(1) 0.304(2) 0.346(2) 0.376(2) 0.391(2) 
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